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Abstract 

We study the dynamical behavior of a single degree of freedom mechanical 
system with a particle damper. The particle (granular) damping was opti- 
mized for the primary system operating condition by using an appropriate 
gap size for a prismatic enclosure. The particles absorb the kinetic energy 
of the vibrating structure and convert it into heat through the inelastic col- 
lisions and friction. This results in a highly nonlinear mechanical system. 
Considering linear signal analysis, state space reconstruction, Poincare sec- 
tions and the determination of maximal Lyapunov exponents, the motion 
of the granular system inside the enclosure is characterized for a wide fre- 
quency range. With the excitation frequency as control parameter, either 
regular and chaotic motion of the granular bed are found and their influence 
on the damping is analyzed. 
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1. Introduction 

In recent years, particle dampers (PD) have been studied extensively 
for use in harsh environments where other types of damping devices are 
not efficient. A PD is an element that increases the structural damping 
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by inserting dissipative grains within holes in a vibrating structure or in a 
box attached to the primary system [1]. The granular particles absorb the 
kinetic energy of the primary system and convert it into heat through inelastic 
collisions and friction between the particles and between the particles and 
the walls of the enclosure. PD are effective over a wide frequency range. 
Moreover, PD are durable, inexpensive, easy to maintain and have a great 
potential for vibration and noise suppression for many applications (see [2] 
and [3]). 

An appropriate treatment of the PD in a given structure requires careful 
analysis and design, parameters such as the size and shape of the particles, 
density, coefficient of restitution, size and shape of the enclosure, and the type 
of excitation of the primary system, among many other factors, are important 
in damping performance (see e.g. |1] and [5]). Due to its complexity, the 
understanding of damping mechanisms involved in PD is still limited. 

The effectiveness of a PD is directly related to cooperative movements of 
the particles inside the cavity. Due to the inelastic collisions and friction, the 
system has a behavior highly nonlinear. The fundamental characteristic that 
affects the performance of the PD is the state of the granular media within 
the PD. Solid-, liquid- and gas-like phases have been identified for vibrating 
granular systems [6]. 

Many authors have focused on studying the behavior of the damping 
device regardless of the complex internal dynamics of the particles. Therefore 
the theoretical models derived from single particle systems [3, IH] are not 
applicable to analyze the nonlinear dynamical multi-particle systems. 

In this work, we use the Discrete Element Method (DEM) [9] to simulate 
the internal dynamics of a PD in order to analyze the nonlinear behavior and 
to characterize the granular system inside the cavity. These particle dynamic 
simulations have been used by others for investigating the behavior of these 
complex granular systems [IHl [HI |12l [13] . The method proposed by Sanchez 
and Pugnaloni [13] to obtain the optimal gap size (space left between the 
granular layer and the top of the enclosure) through the particles effective 
mass is used to achieve the design with the best damping performance for 
the given conditions. 

With the data obtained through simulations, the damping performance 
for a wide range of frequencies is investigated. We distinguish periodic and 
chaotic motion of the granular bed inside the enclosure at different excitation 
frequencies. Results show that the particles run from a periodic motion to 
a chaotic motion. Furthermore, we found that the optimal damping perfor- 
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mance occurs in a window of periodicity. 



2. The mechanical model 

The single degree of freedom (SDoF) primary system (see Fig. [T]) has been 
modeled with a mass M = 2.37 kg, a linear spring K = 21500 Nm~^ and a 
viscous damper with damping constant C = 7.6 Nsm~^. The PD consists of 
= 250 spherical grains in a prismatic enclosure made by six flat walls with 
the same material properties as the particles. The material properties of the 
particles (and walls) and the simulation parameters are listed in Table [l] 
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Figure 1: (a) Model of the SDoF system with a particle damper, (b) Snapshots of the 
particles in the enclosure during a typical simulation. 

The gravitational field g = 9.8 ms~^ is considered in the negative vertical 
direction. Although the SDoF system can only move in the vertical direction, 
the particles move freely inside the enclosure. 

The particles-particle interaction force in the normal direction is based 
on the Hertz-Kuwabara-Kono model [15l [16] . 

-^n = -Ka^^^ - JnVnV<^; (l) 

where = ^E-\J^{1 — v'^)~^ is the normal stiffness (with E the Young's 

modulus, V the Poisson's ratio and = + and Ri the radius of 
particle i), 7n the normal damping coefficient of the viscoelastic contact, and 
i^n the relative normal velocity. 



3 



The particle-particle tangential force Fg, is based on Coulomb's law of 
friction [T5],[l7j. We used a simplified model in which the friction force takes 
the minimum value between the shear damping force and the dynamic friction 

mm- 

Fs = - min {\%VsVa\ , l/Xd^nl) sgn (v^) ; (2) 

where 7s is the shear damping coefficient, Vg the relative tangential velocity 
between the two spheres in contact and fid the dynamic friction coefficient. 
The sign function indicates that the friction force always opposes the direc- 
tion of the relative tangential velocity. 

The system is excited by the harmonic displacement of the base to which 
the spring and viscous damper are attached (see Fig. Let u{t) and z{t) 
be the displacement of the base and the primary mass, respectively. Then, 
the equation of motion for the system is given by: 

Mz{t)+Cz{t)+Kz{t) = Cu{t) + Ku{t) + Fp.^,t{t); u{t) = U cos{ujt) (3) 

where -Fpart(i^) is the z-component of the force resulting from all the interac- 
tions (normal and tangential) of the enclosure walls with the particles. The 
amplitude, U, and the angular frequency, u, of the harmonic vibrating base, 
are control parameters. 

3. Nonlinear dynamic analysis 

The mechanical model described in Section [2] can be numerically solved 
to obtain a scalar sequence of the position along the z-axis of the center of 
mass of the granular bed and the z-position of the box. It has been shown in 
experiments (see e.g. [TH]) that the motion of the granulate is well described 
by the motion of the center of mass, therefore, time series associated with 
the PD are available and these are valid to investigate the granular behavior 
inside the cavity and analyze the damping performance. 

The analysis of nonlinear dynamical systems from time series involves 
linear signal analysis to identify periodic motions, state space reconstruc- 
tion, Poincare sections and the determination of the MLE to identify chaotic 
behavior. We have used for this analysis the Tisean Software package [T9\ . 

The nonlinear analysis requires several inputs for the reconstruction of 
the state space and the determination of the MLE. In this contribution the 
false nearest neighbors method has been used [2D] to determine the minimal 
sufficient embedding dimension and the time delayed mutual information for 
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Shear damping coefficient 


1.098 X lO'' kgs"^m"^/^ 


Excitation amplitude U 


0.0045 m 


Time step bt 


8.75 X 10"* s 


Time of simulation 


196.83 s 


Particle radius 


0.003 m 


Total particle mass 


0.227 kg 


Box lateral side and Ly 


0.03675 111 


Box height 


0.1225 m 


Excitation frequency range 


0.5 to 30.0 Hz 



Table 1: Material properties of the particles and simulation parameters. 



the delay [2T]. The algorithm proposed by Kantz [22] is used to obtain the 
MLE. 



4. Numerical results 



As shown in Table [T| we have simulated the vibration of the system for 
196.83 s. After an initial transient, the system reaches a steady state. A 
simple evaluation of the effective damping yield by the particles can be done 
by fitting the frequency response function (FRF) to a SDoF system including 
only a viscous damper. The amplitude of the response X of a SDoF system 
with no PD is given by 



X = U 



1/2 



(4) 



We carry out a least-squares curve fitting of the DEM data with Eq. |4j 
The values of K, U and C are fixed to the corresponding values in our 
simulations and is the fitting parameter. Figure [2] presents the cavity 
position as a function of excitation frequency compared with the fitting curve 
without PD. 
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Several works [HI [23] predicted a optimal gap value yield the best damp- 
ing performance. We use the particles effective mass to obtain the opti- 
mal gap size [H]. The granular bed at rest has an approximate height of 
Lz = 0.039 m and the gap (L^ — 0.039 m). The FRF is complex for these 
conditions (21]. 




/(Hz) 

Figure 2: Frequency response function of the primary system. The blue continuous hne 
corresponds to the SDoF system without particles in the enclosure and with a mass equiv- 
alent to the system with the PD. The red circles correspond to the FRF for the optimum 
PD. The green diamonds are the MLE for each frequency in the chaotic region. The dark 
gray area corresponds to the transition to chaos, the gray area corresponds to the chaotic 
zone, and the white area to the periodic motion of the trajectory of the center of the mass 
of the granular bed. 

4.I. Periodic behavior 

At low excitation frequency (much lower than the resonant frequency), 
the acceleration of the primary system is less than the acceleration of gravity 
and the granular bed lies always in contact with the floor of the box, thus 
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the center of the mass of the granular bed follows a regular movement. This 
behavior can been seen in Fig. [3j 

In Fig. we plot the trajectory along the z-axis of the center of mass of 
the particles and the motion of the floor and ceiling of the enclosure over a few 
periods of excitation in the steady state regime for an excitation frequency 
/ = 5.5 Hz. We have indicated the position of the granular bed inside 
the enclosure by a band limited by the 2;-coordinates of the uppermost and 
lowermost particle at any given time. To show the regular motion of the 
center of mass at this frequency, we represent in Fig. |3](b) the Poincare map 
(through a stroboscopic map), in this figure can been seen a classical Poincare 
section for a periodic motion of a time series. Moreover, the reconstruction 
of the attractor of the center of mass (see Fig. |8](a)) is typical for regular 
movement, in addition the FFT (Fig. ^c)) shows a peak at the excitation 
frequency of the primary system. 

Clearly, as the response level is too low, there is almost no relative motion 
between the particles and between the particles and the enclosure. Therefore, 
the grains can indeed be viewed as a lumped mass resting on the enclosure 
floor, and the damping effect is very small. In Fig. [2] can been seen that the 
response of the primary system at this low frequency is similar to the system 
without the PD. 

When the excitation frequency increases (/ = 7.5 Hz), the inertia effect 
of the particles exceeds the gravity force and the granular motion develops. 
In Fig. ^d) we plot the motion of the granular bed. This figure shows that 
there is significant particle motion which also leads to impacts between the 
grains and the enclosure floor. The excitation level is not drastic yet and the 
particles do not have enough energy to impact the ceiling of the enclosure, 
then the center of mass continues with a periodic motion. 

The Poincare section (Fig. [sj^e)) is analogous to the one observed for 
/ = 5.5 Hz. The signal has a harmonic spectrum containing only frequency 
components whose frequencies are multiples of the fundamental frequency 
(Fig-gf)). 

4-2. Transition region 

With further increase in / beyond 7.5 Hz a transition phase begins, which 
ends around the frequency / = 11.0 Hz. The granular bed switches between 
quasi-periodic motions and periodic motions. 

For quasi-periodic motions we observe a torus attractor as shown in 
Fig. WCb) for / = 8.0 Hz. The large number of points that lay on a smooth 
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Figure 3: Granular bed response at / = 5.5 Hz and at / — 7.5 Hz. (a) Displacement 
of particles against the enclosure for / = 5.5 Hz. The solid lines show the position of 
floor and ceiling of the enclosure and the colored area indicates the limits of the granular 
bed defined as the position of the uppermost and lowermost particle. The red continuous 
lines is the motion of the center of mass of the grains, (b) Poincare map for / = 5.5 Hz. 
(c) FFT of the trajectory of the center of the mass for / = 5.5 Hz. (d) Displacement of 
particles against the enclosure for / = 7.5 Hz (Colors as in Fig. [3]; a)), (e) Poincare map 
for / = 7.5 Hz. (f) FFT of the trajectory of the center of the mass for / = 7.5 Hz. 
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Figure 4: Granular bed response at f — 8.0 Hz and at / = 10.0 Hz. (a) Displacement of 
particles against the enclosure at / = 8.0 Hz (Colors as in Fig.[3];a)). (b) Poincare map 
for / = 8.0 Hz. (c) FFT of the trajectory of the center of the mass for / = 8.0 Hz. (d) 
Displacement of particles against the enclosure at / = 10.0 Hz (Colors as in Fig. [sjja)). 
(e) Poincare map for / = 10.0 Hz. (f ) FFT of the trajectory of the center of the mass for 
/ = 10.0 Hz. 



9 



closed curve in the Poincare map confirms that the response of the granular 
bed is indeed quasi-periodic, where the response is made up of two frequen- 
cies of which the ratio of one of these frequencies to the other is an irrational 
number. The FFT shows a peak at the forcing frequency, and three other 
secondary peaks; of which one correspond to the particles response frequency 
and other two are simply the sum and difference of the forcing and response 
frequencies. The occurrence of these last peaks, that correspond to irrational 
fractions of the forcing frequency, is evidence that the granular bed response 
is quasi-periodic. 

Alternating between these quasi-periodic motions we observed regular 
motions, which are characterized by clouds of points in the Poincare sections 
and sub-harmonic peaks in the FFT. We show an example in Fig. |4]^e), where 
we find a period-2 motion at / = 10.0 Hz. 

The energy of the particles inside the PD remains low in this region and 
no grains can reach the ceiling of the box. The increase in the amplitude of 
vibration of the primary system produces low granular densities upon each 
impact with the walls of the enclosure and the transfer of the momentum is 
not significant (see Fig. |4]^a) and Fig. |4](d)), therefore the granular damping 
performance is poor and the response of the primary system follows with 
similar characteristic to the system without the PD (see Fig|2]). 

4-3. Chaotic behavior I 

As the value of frequency is increased to 10.5 Hz, wrinkling of the torus 
in the Poincare section is observed. The distortion of the torus becomes 
more significant as / is further increased resulting in the breaking up of the 
continuous curve, forming a fractal torus, at / = 11.0 Hz. Fig. [5] shown the 
behavior at / = 10.5 Hz and at / = 11.0 Hz. 

The response level is high in the range from / = 11.0 Hz to / = 14.0 Hz, 
there is significant relative motion between the particles and between the 
particle and the enclosure with enough energy to impact the ceiling (see e.g. 
Fig. ^g) for / = 13.0 Hz). At these frequencies the Poincare map, which 
consists of a large number of points lying on an open curve (Fig.|5](h)), suggest 
that the attractor of the grains response is strange. This in turn indicates 
that the response could possibly be chaotic. 

The broadband frequency observed in the power spectrum plot (see Fig. [Sj^i)) 
suggests that the motion of particles reaches a chaotic regime. The MLE of 
the grains response for the range between / = 11.0 Hz to / = 14.0 Hz are 
shown in Fig|2j These positive values confirm that the granular bed motion 
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Figure 5: Granular bed response at f — 10.5, at / = 11.0 Hz and at / = 13.0 Hz. (a) 
Displacement of particles against the enclosure at / = 10.5 Hz (Colors as in Fig. [sjja)). 
(b) Poincare map for / = 10.5 Hz. (c) FFT of the trajectory of the center of the mass for 
/ = 10.5 Hz. (d) Displacement of particles against the enclosure at f = 11.0 Hz (Colors 
as in Fig. |3ja)). (e) Poincare map for / — 11.0 Hz. (f) FFT of the trajectory of the 
center of the mass for / = 11.0 Hz. (g) Displacement of particles against the enclosure at 
/ = 13.0 Hz (Colors as in Fig. [sj^a)). (h) Poincare map for / = 13.0 Hz. (i) FFT of the 
trajectory of the center of the mass for / = 13.0 Hz. 
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is indeed chaotic. Moreover, an strange attractor is found for these range of 
frequencies and we show it for / = 13.0 Hz in (Fig. [8|^b)). 



4.4- Optimal damping 

In the range between / = 14.5 Hz to / = 16.0 Hz, we found a window 
of periodicity, where the particles impact the enclosure two times per period 
of excitation. At these frequencies, the granular bed enters a solid-like state, 
colliding both the floor and the ceiling out of phase. We plot this behavior 
in Fig. [6j Note that the compactation of the granular bed is high in this 
regimen and the transfer of momentum at impact also increases. 

The absolute velocities of the center of mass of the particles and the 
primary system are opposite to each other at the instant immediately before 
contact. Due to maximization of the relative velocity, the PD dissipates the 
maximum energy. As shown in Fig|2| the maximum effective damping is 
obtained in this range. We show an example of the periodic attractor at 
/ = 14.5 Hz in (Fig. Mc)). 
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Figure 6: Granular bed response at / ~ 14.5 Hz. (a) Displacement of particles against 
the enclosure (Colors as in Fig. Isja)). (b) Poincare map. (c) FFT of the trajectory of 
the center of the mass. 



4-5. Chaotic behavior II 

Further increase in the control parameter resulted in chaotic vibration in 
the response of the center of mass. This chaotic behavior of the granulate is 
observed from / = 16.5 Hz immediately after the window of periodicity. At 
these frequencies (above the primary system resonant frequency) the motion 
of the grains is very intense and the particles enter a gas-like state. The 
excitation level is drastic and the granular bed expands significantly (very 
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low compactation). Due to the high energy in the system, only a few particles 
hit the walls at each instant and again the systems response is similar to the 
system without the PD (see Fig|2]). 
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Figure 7: Granular bed response at / = 21.0 Hz. (a) Displacement of particles against 
the enclosure (Colors as in Fig. |3|a)). (b) Poincare map. (c) FFT of the trajectory of 
the center of the mass. 

Fig. [7] shows the above resonance chaotic behavior of the trajectory of the 
center of the mass of the granular bed for / = 21.0 Hz. The cloud of points 
in the Poincare section (Fig. ^h)) and the broadband frequency observed in 
the power spectrum plot (Fig. [7](c)) suggest that the response of the grains 
is chaotic in this regimen. Moreover, a strange attractor is clearly identified 
from the reconstruction plot in (Fig. |8](d)). 

Even though the system response presents a chaotic-like characteristic, to 
assure this conclusion we plot in Fig|2]the MLE for this range of frequencies. 

Finally, a comparison of the motion of the granular bed in the chaotic 
regions show that the particles behavior is different for one or other zone. 
The chaotic motion below resonance is characterized by particles colliding 
against the floor and ceiling. For forcing frequencies above resonance, the 
inertia of the granular bed does not allow the center of mass follows the 
motion of the box, thus only a single impact motion with the floor exists. 

5. Summary and conclusions 

In this paper we study a single degree of freedom mechanical system 
with a prismatic PD and we have used the optimal gap size for the best 
performance of granular damping. In particular, we study the dynamical 
behavior of the center of mass of the granular bed. 
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Figure 8: State-space reconstruction, (a) / 
(d) / = 21.0 Hz. 
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We have observed that the grains run from a periodic motion to a chaotic 
motion when the excitation frequency increases. For low excitation frequen- 
cies the particle motion is periodic and the granular damping is very small, 
due to the lack of relative motion between the particles and between the 
particles and the enclosure. 

In contrast, for high excitation frequencies the granular bed has great 
energy and the movement of its center of mass behaves chaotically. Much 
larger frequencies lead the granular bed to expand significantly. The granular 
sample enters a gas-like state with only a few particles colliding with the 
enclosure in each oscillation. This transfers little momentum to the primary 
mass and the system presents low damping performance. 

Is unclear what is the correct transition to chaos, but there are signs 
of transition from quasi-periodicity. The emergence of the chaotic attrac- 
tor in the response is preceded by the alternate appearance of quasi-periodic 
and periodic sub-synchronous vibrations. However, further studies are neces- 
sary to analyze intermediate frequencies to achieve a smoother transition to 
chaos. We found that at frequencies near the resonance of the primary system 
there is a window of periodicity in which the maximum granular damping is 
achieved. This is mainly due to the fact that the particles hit the enclosure 
out of phase and results in a strong reduction of the maximum displacement 
of the primary system. 

We have found that the two chaotic zones have different behavior. Close 
to the resonant frequency, the particles collide against the floor and the 
ceiling. When the excitation frequency increases, the inertial forces do not 
allow the complete motion of the particles. Thus, the granular bed collides 
only against the floor. 
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